Topological Excitations near the Local Critical Point in the Dissipative 2D XY model 
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The dissipative XY model in two spatial dimensions belongs to a new universality class of quantum 
critical phenomena with the remarkable property of the decoupling of the critical fluctuations in 
space and time. We have shown earlier that the quantum critical point is driven by proliferation in 
time of topological configurations that we termed warps. We show here that a warp may be regarded 
as a configuration of a monopoles surrounded symmetrically by anti-monopoles so that the total 
charge of the configuration is zero. Therefore the interaction with other warps is local in space. They 
however interact with other warps at the same spatial point logarithmically in time. As a function 
of dissipation warps unbind leading to a quantum phase transition. The critical fluctuations are 
momentum independent but have power law correlations in time. 



One of the important developments in condensed mat- 
ter physics of the last quarter century is the realization 
that there exists a class of quantum-critical phenomena 
which are not simply extensions [iHj] of the well under- 
stood dynamical critical phenomena near classical critical 
points. In the usual extensions, an essential features of 
classical critical fluctuations that the correlation length 
in the temporal-direction is a power z of the correlation 
length in the spatial direction is maintained. It was ob- 
served by contrast in one part of the phase diagram of 
the normal state of the high temperature cuprate super- 
conductors, that the properties could be understood only 
by quantum-critical fluctuations which are local spatially 
but with a simple (oc 1/t) power law correlation in time 

Subsequently similar local criticality has been ad- 
duced from measurements of the correlation functions 
for some heavy-fermion compounds near their antifcrro- 
magnetic quantum-critical point It has also been 

sug gest ed that disordered superconducting thin films 



[lCt and superconducting wires [12J, [13| have local 
quantum-criticality at the superconducting to insulating 
transition [I3 - [l7j . Local phase slips have been shown to 
play a crucial role in dissipation driven transitions in one 
dimensional superconducting wires [l8| . The wire can be 
thought of as a linear chain of coupled superconducting 
grains. Each event in space-time creates a dipole consist- 
ing of a phase slip in one junction and an anti-phase slip 
in a neighboring junction. 

In all the physical situations mentioned above, one may 
identify regimes in parameters like temperature, doping, 
pressure, thickness and magnetic field or charge density 
where the global properties are dominated by local de- 
grees of freedom. The absorptive part of the fluctua- 
tion spectrum of these degrees of freedom displays uo/T 
scaling. This is simply the thermally weighted Fourier 
transform of 1/t correlations. There is no evidence of a 
spatial correlation length which is related to such a tem- 
poral correlation. 

For the dissipative 2DXY model, a variational and 
weak-coupling calculations suggested that the quantum- 
criticality is local 14 , [HI . To compute the correlation 



functions, a non-perturbative treatment is necessary. Re- 
cently, such a theory of this and some related models has 
been developed 19,[20||. The principal accomplishment 
is to show that spatial and temporal degrees of freedom 
of fluctuations exactly decouple in the partition function 
because the fluctuations can be represented in terms of 
orthogonal topological defects - vortices which are shown 
to interact locally in time and logarithmically in space 
and a new class of defects warps, which interact locally 
in space and logarithmically in time. 

In this paper we resolve some remaining questions 
about the discrete nature of the warps and explain why 
the warps interact locally in space. We do this by giving 
an explicit representation for them in terms of config- 
uration of phase slips for the 2d dissipative xy model 
and give a physical picture of how they are orthogonal 
to the vortices. This is the only two-dimensional model 
for which local quantum-criticality and uj/T scaling have 
been explicitly proven. As such, the general understand- 
ing and technical procedure may be important for other 
problems where similar criticality appears to govern. 

This paper is organized as follows. In the next section, 
we summarize briefly the transformations 3, 13] leading 
to the Action which displays decoupling of space and 
time variables. In the subsequent sections, we introduce 
Vortices and Warps, give their physical description in 
terms of configurations of phase slips, and show that the 
Action in terms of them is asymptotically identical to the 
Action derived earlier. 



DISSIPATIVE 2DXY MODEL AND ITS DUAL 
REPRESENTATION 



The classical 2d xy model consists of U(l) degrees of 
freedom, represented by an angle 9, living on the sites (ij) 
of a regular lattice, assumed here to be a square lattice 
see Fig. (JlJ, with a nearest neighbor interaction of the 
form 
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FIG. 1: The directed link variables are labelled as shown in 
the (a), (b )We define a two component vector living on the 
sites of the original lattice whose components are the two 
directed links variables: m = (mij-i +x j,rrii t j-i t j +x ) 



Since a continuous symmetry cannot be spontaneously 
broken in two dimensions [U [22| , this model does not 
support a long range ordered phase. Nevertheless a phase 
transition does occur at finite temperature where the 
correlation function of the order parameter e changes 
from exponential to power law. This is the Kostcrlitz- 
Thouless-Berezinskii transition [2^, [24]]. The quantum 
dissipative generalization of the model includes two dy- 
namical terms and is given by 



Z = f DBi (r) exp - j* dr I £ £ {d T 6%) - J £ cos 
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where C is the capacitance and a = Rq/R where Rq = 
h/Ae 2 . 

The physics of this phase transition is better under- 
stood in terms of the topological defects of the system. 



To do so we follow the standard procedure of using the 
Villain transform and integrating out the phase degrees 
of freedom [20] . The Villain transform involves expand- 
ing the periodic function in terms of a periodic Gaussian 



exp 
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where rriij : ki are integers that live on the links of the original square lattice. We can combine the two link 
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variables mij-i+ij and TOij ; jj+i into one two component 
vector mi j that lives on the site {i, j} of the lattice (see 
figH]). We expand the quadratic term and transform to 
Fourier space. Keeping the leading quadratic term Q± j — 
~ —a^x^xy, where a is the lattice constant, x — ai 
and y — aj, we get 

(0 xy - 9x+i,y - 2vrm* J w a 2 V 2 x xy + AnaV x 6 xy m Xty 

+ 4^ 2 m^, (4) 

where m% y is the x component of the vector field given 
by the integer m X y X +i y . In the absence of dissipation, 
there is competition between the kinetic energy and the 
potential energy terms, the former minimized by a state 



where fty is disordered stabilizing an insulating phase 
while the latter minimized by a fixed value of Oij stabi- 
lizing a superconductor. Since 0's are Bosonic degrees 
of freedom, we need to impose the boundary condition 
that Oij (/?) = 9ij(0). The periodicity in the imaginary 
time direction and the compactness of the field Oij im- 
plies that there is an additional degree of freedom that 
has to be accounted for which is the winding number. At 
T = 0, the (imaginary) time direction becomes infinite 
in extent, and the non-dissipative model is in the 3DXY 
universality class. First we discretize the imaginary time 
direction in units of At and work on a three dimensional 
lattice. Introducing the variables m which only live on 
the spatial links, the action is [3, l20| 
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where c = a/ At, J — > Ja 2 Ar, C — > Ca 2 /Ar and 
m — > raj a. We have also redefined a — > aa 3 . The 
last term in the denominator is unimportant and may 
be dropped. The action in Eqn(S]has two possible phase 
transitions, depending on whether the capacitance C or 
the dissipation term oc a in the numerator of the second 
term in Eq. ([5]) dominates in long wavelength and low 
frequency limit. The former corresponds a critical point 
with the dynamic critical z = 1, i.e. we recover the loop 
gas model which belongs to the 3d xy universality class. 
The latter corresponds, as we will show to z — oo, i.e. a 
fixed point with local criticality. The two terms have the 
same scaling form for z = 2. 

Suppose we separate m into the usual transverse part 
in terms of a vortex field p v (k,u n ) 



ik x m t (k,w„) = p v (k,w n ) 



(6) 



and a longitudinal part through introducing the warps 



w„m£(k,u;„) = ckp w (k,0J n ). 



(7) 



Quite miraculously the partition function in Eq. ([5]) can 
be written exactly as 



= ex p 
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where 



G(k.u n ) 



Jck 2 + (C/c) ul + a |cj„ | k 2 



(9) 



The first terms describes vortices interacting logarithmi- 
cally in space with no retardation in time, the second 
describes warps interacting logarithmically in time with 
local interactions and the last term is irrelevant near the 
critical point. The partition function near criticality can 
then be easily be evaluated. This is the procedure fol- 
lowed in Ref. sua to evaluate the correlation func- 
tions near criticality. 

Two questions immediately arise. (1) Is the discrete- 
ness of the variables mi being respected in this proce- 
dure? While this is not obvious even for a vortex in the 
continuum equation for vortices, Eq. ([6]), it is well under- 
stood that such equations indeed respect the discreteness 
of m t through the discreteness of p v . A less than satis- 
factory answer for the discreteness of the variables 
was provided in Ref 20. A clear demonstration is pre- 
sented below. (2) While the locality of the interaction 
between warps is explicitly shown in Eq. (JHJ) , it is worth- 
while to show this in direct physical terms. This is also 
accomplished below. 
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Vortex 

FIG. 2: An isolated vortex at site fi. In terms of the vector 
field m, the vortex is a string (dashed line) that stretches out 
to the system boundary. All links intersected by the string 
have a constant m%j represented by the red bonds. 



TOPOLOGICAL DEFECTS 
Vortices 

The 2Dxy model supports topological excitations 
where the vector field m acquires nonzero curl. An iso- 
lated vortex, of unit vorticity, at site y, is represented in 
the figlU In terms of the dual variables m, the vortex is 



a nonlocal object. A string stretching from site /x (greek 
alphabets are used for sites of the dual lattice) to infinity 
is a representation of the vortex with the link variables 
being finite only on those links intersected by the string. 
If {xirp,,yirfj,) labels the lower right hand corner of the 
plaquette labeled /i, we can represent the vortex as 



m v (x,y)=yQ(x — xi r , i )6(y-yi ri j,) (10) 

where a is the lattice constant. Note that V x m = 
S(r— Ti) . An ensemble of such strings faithfully generates 
all possible configurations of vortices. The vortex degrees 
of freedom are defined as V x m (r) = p v (r). 
The classical partition function is completely determined 
by the vortex degrees of freedom. The 2DXY model is 
mapped to a neutral plasma of vortices (zero total vortic- 
ity). As temperature is raised through the critical tem- 
perature, vortices unbind and proliferate. As a result the 
order parameter correlation function changes from power 
law at low temperatures to exponential at high temper- 
atures. 

Dynamics: The quantum model includes the dynamics 
of the vector field m. The time evolution of the vector 
field results in a vortex current given by 



J v = z x — (11) 

(XT 

The vortex density and current satisfy the continuity 
equation. In the absence of dissipation, the partition 
function is 



Z = ^]exp 
m 



1t 2j Jc \Pv ( k ^«)| 



(C/c) \3 V (k,w„)| 



(C/c) ojI + Jck 2 



(12) 



Eq. (|12p shows the well known fact that the quantum 
2DXY model without dissipation maps onto a three di- 
mensional vector field with coulomb interaction and a 
constraint of no divergence (continuity equation). In ad- 
dition to spatial phase slips on links this model allows for 
phase slip events in time. 

Since the current involves a change in the value of the 
vector field, consider an event in time where a single 
bond acquires a nonzero value. Such an event produces a 
vortex and and antivortex on neighboring plaquettes as 
shown in fig|3] 

Alternatively, one can interpret the event of creation of a 
vortex dipole as an vortex current. In time interval St a 



vortex moves from site 1 to site 2 in figUl leaving behind 
an antivortex at site 1. The total vorticity before and 
after the event is zero as is required by the condition of 
having a neutral plasma. Over time St a phase slip event 
occurs on the link between sites 1 and 2 at time Tj. If 
the lower end of the bond is labeled by v xy = {x, y}, the 
event, in our notation, is given as 

m = -yS (r - r xy ) 6 (r - n) (13) 
The vortex current generated by such an event is 

Jy =xS(r- r xy ) 5 (r - n) (14) 
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Phase Slip event on a link 

FIG. 3: A phase slip between nearest neighbors in time pro- 
duces a vortex and an antivortex on neighbouring plaqiettes. 
In time one bond has acquired a nonzero value, represented 
by the red bond. If the value on the bond is -1, the vorticity 
at site 1 is negative and on site 2 is positive. 




FIG. 4: A phase slip event in rime results in a change in the 
link variables. For a change of 2tt at site {i, j}, the four links 
connected to it acquire the values shown in the figure. 



Thus a phase slip on a link is equivalent to a local in 
space and time vortex current. 

Local phase slips and Warps 

The periodic in time boundary condition allows for 
phase slip events on a site. Such events change the 
winding sector and dynamics of the vector field that is 
not captured by the vortex current. Consider the ef- 
fect of a change of 2tt at site {i,j} so that 0y winds 
around to #y 4- 2tt over a time St. Then the four spatial 
links connected to the site experience phase slips. 

As shown in Fig. ((4]), the corresponding link variables 
are: m hJ . l+xA = 1, m JJ:JJ+1 = 1, ///,.,., i . ; = -1 and 

^i,j:i,j—X 1- 

The phase slip event at time r, and site has the fol- 
lowing representation in terms of the vector field m: 

m(r,r) = [($ + y) S (r - ry ) - xS (r - r# - osc)(15) 



- yS(r- Vij - ay)] Q (r - n) 

Such a vector field distribution has no curl and hence does 
not effect the vorticity. On the other hand the divergence 
is nonzero and phase slip events generate field configu- 
rations that are orthogonal to those created by vortices. 
For a general vector field one does expect two kinds of 
sources to generate an arbitrary distribution. For the 
2+1 dimensional quantum model we have, besides the 
vortices, the additional topological entity to describe the 
winding number sector in time. Events that change the 
winding number sector, i.e. local phase slips, acts as 
sources for a divergence in the vector field. Just as a 
vortex is equivalent to an electric charge in the dual lan- 
guage, the sources created by phase slips can be shown 
to be a local distribution of monopoles (p m )- Given the 
distribution in eqn ll5[ the corresponding configuration of 
monopoles, which we term the charge of the phaseslip is 

p m (r,r) = V-m(r,r) (16) 

= [4S (r - v^) -8{v- r y + ax) 

- 5 (r - v i: j - ax) - 8 (r - v l} + ay) 

- S(v-Vij -ay)]Q(r-Ti) 

The monopole distribution equivalent to a phaseslip is 
shown in figO The total monopole charge of the con- 
figuration is zero. Since the distribution has azimuthal 
symmetry all harmonics are zero. This is the two dimen- 
sional lattice realization of the configuration of a charge 
surrounded by an equal but opposite charge distributed 
over a spherical shell of radius a in three dimensions. The 
magnetic field due to the charges is confined within one 
unit cell around the site of the phase slip and is zero out- 
side. Thus two phase slip events can interact only if they 
are at most one lattice spacing apart; the interaction is 
local in space. Although this is physically obvious from 
this discussion, we will demonstrate this explicitly in the 
next section. 

Phase slip events generate a local vortex current which 
is divergenceless but has a finite curl. 

J ps (r, t) = [(-x + y)8(v- r^) - yS (r - r y - ax) 

+ x6(v- v^ - ay)} S (r - n) (17) 

The continuity equation is satisfied as this current has 
no divergence. 

DISSIPATION DRIVEN TRANSITION 

Any configurations of the vector field m can be de- 
scribed by an ensemble of phase slips and vortices. Since 
the former are divergence free and the latter curl free, we 
have 
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FIG. 5: Warps are configurations of the tieid m 
divergence but no curl. For a warp of unit strength at site 
{i, j}, the divergence has a magnitude of 4 at the site and -1 
at the four nearest neighbor sites. In general the magnitude of 
the divergence at the site will equal to the number of nearest 
neighbors on the lattice. 



Vxm t (r,r) = p v {r,r) (18) 
m/(r,r) = ^ p l ps G (t - n) [(x + y) 8 (r - r y ) 

i 

— xS (r — Yij ~ ax) 

- yS(r~ - ay)] 

I 



The locality in space is reflected in eqnQ~9] by the 
factors (1 — exp (— ik x a)) and (1 — exp (— ik v a)). Since 
k x a,k y a <C 1, to leading order we get 

m, (k, W „) « 1^M±Ml p (k)Wn ) (20) 

The two components of the field m t and m; are orthog- 
onal to each other. Comparing the longitudinal compo- 
nent of nij generated by a phase slip in the long wave- 
length limit and a warp (eqnJT]), we find they are pro- 
portional to each other. In particular kap ps — p w . The 
partition function can be recast in terms of the charges. 



Z = ^ exp 

{Pv,P V s} 



E -4^12 \P« (k^n)| 2 -47T 2 



(4 - 2cos(fc x a) - 2cos{k y a)) \p ps (k,w„)| (21) 



c\uj n \ 



4tt 



JC/c-a\u n \C/(?-a 2 k 2 /c 
(C/c) u)% + Jck 2 + a \ui n \ k 2 



(4 - 2cos(k x a) - 2cos{k y a)) \p ps (k,w„)|' 



The singularities decouple space and time. Vortices in- 
teract logarithmically in space with no retardation while 
phase slips interact logarithmically in time. The interac- 
tion between phase slips is short range in that two phase 
slips interact only if they are less than a lattice spac- 
ing apart. As pointed out in the previous section, phase 
slips are a localized configuration of monopoles whose 
magnetic field is confined to one unit cell. Thus the in- 
teraction is short ranged. 

The form of the action is similar to Eq. JS]), which was 
the basis for the calculations of Ref . [19l |20( ■ The only 
difference is in the factor of (4 — 2 cos (k x a) — 2 cos (k y a)) 
that multiplies the terms involving warps, which comes 
from explicitly deriving the action on a lattice as opposed 
to in a continuum in Eq. ([5]). The nature of singularity 
is not altered. To make this point explicitly, consider the 
long wavelength limit of the most singular part of the 
action for the phase slips (second term in eqn l2~T|) . For 



k x a, k y a <C 1, this term is (upto factors of 47r 2 ) 

— r (4 -2 cos (k x a) - 2 cos (k y a)) \p ps (k,w n ) 

c\u> n \ 

~ , a , k 2 a 2 \p ps (k,<^„)| 2 = " \p w (k,uj n )\ 2 (22) 
c\ui„\ C\0J n \ 

The fact that the discrete phase slips interact locally in 
space translates to a local interaction among warps in 
the long wavelength limit. Thus the critical theory are 
identical. At the microscopic level, the connection with 
phase slips guarantees that the discrete nature of the m 
fields is obeyed. 

The last term in eqn l21l is not singular as it involves a 
scalar field with three dimensional Coulomb interaction. 
The most singular interaction among the phase slips in 
real space can be written as 
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ex p 

{pp.} 



y^J dr J dr'Aapps (R 4 , r) log \w c (r - r')| ,0p S (Rj, r') 
^ J dr J d,Tap ps (R 4 ,r)log|o; c (t - t')\ p ps (R 4 + r 3 ,r') 



(23) 



where r.,- runs over the four nearest neighbors of i on 
the square lattice and oj c is the high frequency cutoff. 
Furthermore the singularity at uj — implies that at 
each site the sum of all charges over time is zero. In 
addition to on site long range in time interaction, phase 
slips involve interaction among nearest neighbors. As a 
there is a unbinding transition in time at a c = 1/4. The 
phase transition is identical to the one analyzed in ref.[ 
19L 2o| . The correlation functions are local in space and 
power law in time. 



CONCLUSION 

The local character of the quantum phase transition of 
the dissipative 2DXY model is a consequence of a new 
topological defect, warps, that arises due to phase slip 
events in time. The discrete nature of the charges is 
related to the periodic in time boundary condition en- 
forced on bosonic fields in the path integral formalism. 
The compactness of the phase field 9 allow for world lines 
in time that possess nontrivial winding numbers. The 
events in time where the winding numbers change cor- 
respond to spatial configuration of phase slips between 
nearest neighbor links that are orthogonal to the vortices. 
While in dissipation free systems such charges contribute 
nonsingular terms to the effective action, ohmic dissipa- 
tion yields a new class of quantum critical point. 

The universality class defined by such singularities is 
a new direction in the study of quantum critical phe- 
nomena with potential realizations in a number of sys- 
tems such as cuprate superconductors, disordered two 
dimensional thin films, heavy fermion superconductors 
and possibly the pncitide superconductors near their an- 
tiferromagnetic quantum-critical point, as well as plateau 
transitions in quantum Hall effects. We speculate that 
these systems also belong to the universality class where 
the spatial and temporal correlations are decoupled due 
to orthogonal topological defects. The precise nature of 
the defects is bound to have different microscopic origin 
and should be the subject for future studies. 
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